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We analyze the initial value problem for scalar fields obeying the Klein-Gordon 
equation. The standard Cauchy initial value problem for second order differential 
equation is to construct a solution function in a neighborhood of space and time 
form values of the function and its time derivative on a selected initial value sur- 
face. On the characteristic surfaces the time derivative of the solution function 
may be discontinuous, so the standard Cauchy construction breaks down. For the 
Klein-Gordon equation the characteristic surfaces are null surfaces. An alternative 
version of the initial data needed differs from that of the standard Cauchy prob- 
lem, and in the case we discuss here the values of the function on an intersecting 
pair of null surfaces comprise the necessary initial value data. We also present an 
expression for the construction of a solution from null surface data; two analogues 
of the quantum mechanical Hamiltonian operator determine the evolution of the 
system. 
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1. Introduction 

The standard Cauchy problem is to use a partial differential equation and the 
values of the solution function and some of its derivatives given on an initial value 
surface to determine the values of the function in a neighborhood of the surface 
[1] The derivatives needed for an n th order partial differential equation in time 
are the first n — 1 time derivatives, and in applications to physics, a surface often 
used is three-space at some initial time. Solution of the standard Cauchy problem 
proceeds by showing how initial value data and the differential equation determine 
all the time derivatives of the function on the initial surface and thus allow a power 
series development of the function for future times. 

If the solution function has a discontinuity in its n — 1 time derivative across 
some special surface the the standard Cauchy problem cannot be solved using that 
surface for the initial value data. Such special surfaces, the characteristic surfaces, 
or simply characteristics, have been studied extensively for many of the differential 
equations of physics. The physical significance of characteristics is that the solution 
or a derivative on one side may differ discontinuously from that on the other side, 
so the characteristic can represent a physical discontinuity interpretable as a wave 
front or shock wave. For example, the characteristics of the Maxwell equations in 
free space are the wave fronts of light, that is null surfaces [2]. The characteristics 
of the Einstein equations in free space are also null surfaces, so that gravity may 
be interpreted to propagate at the velocity of light [2], [3], [4]. 

The characteristics of the Klein-Gordon equation for free particles in flat space 
are easily shown to be null surfaces, as are the characteristics of a large class of 
interacting versions of the Klein-Gordon equation. Thus the derivative disconti- 
nuities of the solution functions propagate at the speed of light, independently of 

3 



the mass parameter in the equations. Note in particular that the discontinuities 
do not propagate at either the phase or group velocity of the wave solutions. 

For many years, the use of null coordinates - coordinates lined up with the null 
surfaces - has been common in general relativity. For example, many problems 
involving black holes are conveniently studied with null coordinates since black 
hole surfaces are time independent null surfaces [6], [7]. Recently interest has 
developed in using null coordinates in the study of quantum field theories such as 
quantum electrodynamics and quantum chromodynamics this is often referred to 
as light-cone quantization, although the planar surfaces x — ±i are usually used 
[8]. In quantum field theories the behavior of the field operators is related to the 
behavior of the classical (c-number) fields. For example, the propagator may be 
interpreted as either the Green's function of the classical field or as the vacuum 
expectation value of the time-ordered field operators [10]. It is thus possible that 
the discussion of this paper may have implications for studies of black hole physics 
in curved space or quantum field theories in flat space. Regarding the latter topic, 
the transition from classical Green's functions to quantum vacuum expectation 
values is beset with subtleties stemming from the possible appearance of anomalies 
upon quantization [9]. 

We have studied the Klein-Gordon equation and its characteristics, in particu- 
lar how the values of the solution function on characteristics determine the values 
of the fields in a neighborhood of space and time. Operators that determine how 
the solutions develop in time are usually called evolution operators, and we use 
here the same term to include similar operators associated with null coordinates. 

We first briefly review the standard Cauchy problem for the Klein-Gordon 
equation, and give an operator expression for developing the solution from a con- 
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stant time surface using as initial value data the function and its first time deriva- 
tive. We then obtain the characteristics and solve the analogous problem for the 
planar null surfaces x = ±t. Initial data for this situation consists of the value of 
the function on the pair of null surfaces, and does not include derivatives. We 
give a simple evolution expression in terms of null coordinates, which is a direct 
analogue to the usual expression of quantum theory, except that two analogues 
of the Hamiltonian operator determine the evolution solution [11], [12], [13], [14], 



The standard form of the Cauchy problem for the Klein Gordon equation is 
to determine how a solution may be developed from the appropriate initial value 
data on a constant surface, which we take to be t = 0. (See figure 1.) The Klein 
Gordon equation for the scalar function <p(x, t) is 



(We use the units in which h = c = 1. The slash notation indicates differentiation 
with respect to the indicated variable, in this case time t.) 

To solve the Cauchy problem we use the Klein Gordon equation and the value 
of the solution function and its first time derivative on the surface t — to compute 
the value of the solution at future times. For this we use a Taylor series expansion 
in t, which we write as 



[15]. 



2. Development of a Scalar Function 



From a Constant Time Surface 




(2.1) 




(2.2) 
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We specify the solution and its first time derivative at t = as 



<f)(x , t) t=0 = h(x) <j)\ t (x,t)t=o = k(x) (2.3) 

Then the Klein Gordon equation allows us to calculate all the remaining time 
derivatives on the t — surface needed in the series expansion; the second derivative 
is 

<f>\t\t& t) t =o = [V 2 </>0?, t) - m 2 (j)(x, t)]t=o = [V 2 - m 2 ]h(x) (2.4) 
Higher derivatives are obtained in the same way, as 



<t>\2nt{x,t)t=o = [V 2 - m 2 ] n h(x) = (-l) n E 2n h(x) (even or 2n derivative) 



0|(2n+i)t(^*)*=o = [V 2 -^ 2 ] n ^(^) = {-l) n E 2n k(x) {odd or 2n+l derivative) 

(2.5) 

Here E 2 is the Klein Gordon energy squared operator, E 2 = m 2 — y 2 . Now 
the series (2.2)may be written as 



This may be formally summed to yield a concise result, 
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<j)(x, t) = cos(Et)h(x) + ^k(x) (2.7) 

E 

which displays the solution in terms of evolution operators acting on the initial 
value functions h(x) and k(x) given on the initial t — surface. 

The formal result (2.7)is directly usable if h(x) and k(x) are expressed as 
Fourier transforms; the operator E may then be replaced by the value ±^/m 2 + p 2 , 
for the p-th Fourier component. For the p-th component the appropriate h(x) 
is e tp - x , and the appropriate k(x) is —iEe ip - x . Substituting these into (2.7)one 
obtains the correct function e - lEt + l P-x f or foe p-th component. It is thus clear 
that (2.7)indeed represents the general solution of the Klein Gordon equation. 

So far we have limited our discussion to the free Klein Gordon equation (2.1). 
However if there is an additional interaction term in the equation which is a function 
of the scalar field and time and position, F(<f>, x, t), then it is evident that the same 
considerations allow one to determine higher time derivatives just as in (2. 4) and 
(2.5), so that a power series solution will also exist in this more general case. 

3. Characteristics of the Klein Gordon Equation 

Surfaces other than t — may be used for the series development of the Klein 
Gordon solutions; indeed only a limited class of surfaces does not allow such a 
series development - the characteristic surfaces. On these surfaces the standard 
Cauchy problem breaks down. 

Let us attempt to develop a solution of the Klein Gordon equation from a 
surface S, given by t — T(x). (See figure 1.) We suppose that the solution and its 
first time derivative are given on S by initial value functions 
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<f>{x, T(x)) = <j)(x, t) t =T = h(x), (/>\ t (x, t) t=T = k(x) (3.1) 

For brevity we write t = T to denote t = T(x). We then represent the solution 
as a Taylor series in the form 



(x, t) = cf>(x, t)t= T + [t- T{x)]4> lt (x, t)t= T + [t T } X)] <t>\t\t& t)t=T + - (3.2) 



A critical step is to find the second time derivative in this expansion. To do 
this we first relate derivatives of the initial data functions to derivatives of the 
solution function (ft on S. From (3.1)we obtain immediately 

\/h(x) = V0(£, t)t= T + <t>\ t (x, t) t =T V T(x) 

V k(x) = V 0|t(x , t)t= T + <l>\t\t(x, t)t= T V T(x) (3.3) 

where the gradient on <p represents the partial with respect to the 3 spatial 
variables. We differentiate this again and rearrange terms using (3.3)to obtain 



V 2 h(x) = V 2 0(x, t) t=T +(f> lt (x, t) t = T \j 2 T(x)-(t) W {x, t) t=T s 7 T(x) 2 +2\/k(x).S7T(x) 

(3.4) 

Using (3.4), the Klein Gordon equation, and the definitions (3.1)for the initial 
value functions we get the desired equation for the second time derivative of on 

S, 
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(l- V T(f) 2 )0|,| t (f,t)t = r = (v 2 -rn 2 )h(x)-k(x) V 2 T(x)-2 V k(x). V T(x) (3.5) 



If the square of the gradient of T is not equal to 1 we may solve this as in 
section 2, and develop a power series solution in a similar manner. Indeed the 
previous case corresponds to the function T having zero gradient and Laplacian. 
However If the gradient squared is equal to 1 then the equation cannot be solved 
for the second time derivative and the surface is a characteristic. Thus the equation 
for the characteristics of the Klein Gordon equation is 

V T(f) 2 = 1 (3.6) 

It is easy to see that this is the equation of a null surface - that is the normal 
to S is a null vector; to do this we write the equation of S as 

t - T(x) = G(x a ) (3.7) 
The normal to S is the 4-gradient of G, or 

n a = G la = (l,- V T(x)) (3.8) 

Thus n is a null vector by (3.6), and S is a null surface (figure 2). 

In the next sections we will study the particular null surfaces for T taken as 
±x, or in terms of the commonly used null coordinates u and v (figure 3), 

9 



u = t — x = u characteristic or null surface 



v = t + x = v characteristic or null surface 



(3.9) 



It is interesting that the characteristics are independent of the mass m. This 
is an example of a more general fact, that if there is an interaction term of the 
form F(4>, x, t) in the Klein Gordon equation then an equation analogous to (3.5)is 
obtained, with the same left hand side. Thus the equation for the characteristics 
is independent of the extra term and the characteristics are again null surfaces. 

Equation (3.5)gives the equation of the characteristics, but it also tells us that 
the functions h and k on the characteristics are constrained if the solution is regular, 
since the right side of (3.5)is zero. Such constraints are typical of characteristics, 
and in this case there is a simple physical interpretation. To see this consider the 
characteristic t = x, and let be a plane wave of the form e ~ lEt + kx . Then setting 
the right side of (3.5)to zero we see that E 2 = k 2 +m 2 , the usual energy-momentum 
or mass shell relation. 

4. Construction of a Solution Function From Null Surface Data 

We now show how to construct solutions to the Klein Gordon using initial data 
on null surfaces; the pair of null surfaces u = t — x = and v = x + t = will be 
used. In terms of the null coordinates u and v the Klein Gordon equation is 



Here the perpendicular symbol y_|_ refers to the y, z subspace; in the remainder 




(4.1) 
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of this paper we will suppress that subspace dependence and write the Klein Gordon 
equation in the form 



(P\u\v(u,v) = -fl 2 (p(u,v), V 2 = — ( 4 - 2 ) 

For the case of fx = this is the wave equation and has particularly simple 
solutions - any differentiable function of u or v. We accordingly expand about 
this special massless case by setting up a solution as a power series in fx 2 , written 
explicitly as 

(f>(u, v) = (f)(°\u, v) + fi 2 (p {1) {u, v) + fx 4 (/) {2 \u, v) + ... (4.3) 

We substitute this into (4.2)and equate powers of fx 2 to obtain a perturbative 
solution. The zeroth order equation and solution are 

<^ = ' <t> {0) (u,v) = f(u)+g(v) (4.4) 
where / and g are differentiable functions. The first order equation is 

*S, = -* ( °W) = "[/(«) +9(v)] (4.5) 

This is an inhomogeneous linear equation, so the solution is composed of a 
homogeneous part and a particular part. But the homogeneous part is of the same 
form as the solution of the zeroth order equation (4.4), and can be absorbed into 
that solution . Thus we need only obtain a particular solution to (4.5), which is 
elementary. We write it in the form 
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U V 

<jP\u,v) = -[J f(u')du'](v-v )-[J g(v')dv'](u-u ) (4.6) 

Uo Vo 

The parameters uq and vq are arbitrary constants and correspond to arbitrari- 
ness in the functions / and g. 

Higher orders are handled in the same way. The second order solution is 



WW V V 

^ 2 \u,v) = -[Jdu'J f{u")du"\^^ -[jdv> J g(v'')dv'']^-pl (4.7) 

WO Uo Vo Vo 

The series solution can be written as 



' n\ z — ' ml 

n=0 m=0 

where the multiple integral operator r™ on f is defined as 



w u 



Tlf{u) = J du' J du"... J du n f{u n ) (4.9) 

WO W() Wo 

The operators T^ v are linear and the series (4.8)may be summed to give a 
remarkably concise result in terms of a pair of exponential evolution operators, 



v ) = e-^-^/H + e -^ {u - u ^ T ^g{v) (4.10) 

This shows explicitly how the solution is given in terms of the two generating 
functions f{u) and g(v), right moving and left moving waves which are related to 
initial value data as we will discuss below. 
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We illustrate the simplest case of using (4.10)and consider the evolution op- 
erator acting on a constant f(u) = 1. The terms in the series (4. 8) are trivially 
obtained by integration and give 



= E - LJ =Fc(uv) (4.11) 

where we have taken uo = vo = for convenience without loss of generality. 
This is easily verified to be a solution of the Klein Gordon equation by substitution; 
moreover we recognize the series as the Bessel function 

F c {uv) = Jo^v 7 ^) (4.12) 

This special solution will be useful in the discussion of exponential solutions 
below. 

Consider the values of on the null surfaces u = uq = and v = vq = 0. We 
see that 

0oM) = /(«) + <7(0), <hM=f(0)+g(v) (4.13) 

Thus up to a constant the generating functions are the initial values on the 
pair of null surfaces, representing right and left moving plane waves. The value of 
either of the generating functions at may be chosen arbitrarily, and we choose 
for convenience both equal to ^ at the intersection of the surfaces, so that 
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f(u) = 0o(u, 0) — , g(v) = O (O, v) — (4.14) 



We emphasize that the initial value data needed for the null surface develop- 
ment consists of the function on a pair of null surfaces; it is not the same as the 
data needed in the standard Cauchy problem, which is the function and its first 
time derivative on an initial surface [14], [15]. 



5. Plane Wave Solutions 

We next discuss plane wave solutions of the Klein Gordon equation, which 
form a convenient set for the expansion of a general solution. We will display 
the appropriate generating functions and verify that the evolution operators in 
(4.10)give a correct solution. Plane wave solutions of energy E and momentum k 
may be written in terms of t and x or u and v coordinates as 



<t>{t, x) = e ~ iEt+ikx = e -^u-trt; ( 51 ) 

where the null momenta A and r are related to E and k and the mass parameter 
H 2 by 



. E + k E-k x E 2 -k 2 m 2 2 

A = — 7. — , T = — ^ — ,Ar= = — = ur (5.2) 

2 2 4 4 



For the plane waves the appropriate generating functions are thus 
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f(u) = e-* Xu - 1 -, g(v) = e- i ™ -± At = f (5.3) 

We substitute these generating functions into (4.10), abbreviating for conve- 
nience B = —iX, C = —it, jj? = —BC, and find 



<t>(u,v) = (e BCvT "e Bu + e BCuT "e Cv ) - (e BCvT " + e BCuT «)- (5.4) 

The last term in this we recognize as F c which we obtained as the series in 
(4.12). The first two terms are readily evaluated by expanding in double series as 



BCvT Uo Bu^ BCuT v „Cv,_ ^ (BCvT u ) n (Buy | ~ (BCuT v ) m {Cv) k 

n\ j\ ^ m\ k\ 



/ BCvF Up Bu + p BCuT Vp Cv\ yu^ui. u j ^ua r 

n=0,j=0 J m=0,fc=0 



(5.5) 



It is easy to obtain the operation of T n on powers of u, as 



ui +n 



Ku J = T. 7TT- 7 (5-6) 

(j + l)(j + 2)...(j+n) 1 ; 

We substitute this into (5.5)and rearrange summation indices to find 



(p BCW Up Bu , p BCuT Vp Cv, _ ^ {Cv) n W A (Cv) m {Buf 

1 + j 2^ n \ j\ + ^ ml k\ 

n<j m>k 



V v " 7 v , r + E lz ri^ - e Bu+Cv + fiH (5.7) 



(C< y, (V^) 2 _ „Bu+Cv 

n\ ?! ^ 

n=0,j=0 J n=0 

Thus the F c part cancels and in summary of (5.4)to (5.7)we have 
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(f)(u, v) = e Bu+Cv = e ~ lXu ~ iTV (5.8) 

which is the correct plane wave solution. It follows that the evolution operator 
equation (4. 10) will produce the correct solution in any case that may be expanded 
as a Fourier transform. 



6. Relation to Null Surface Quantization 

The form of the solution in (4.10)is suggestive of the quantum mechanical time 
evolution operator e~ %m . In fact the analogy can be made complete and is the 
basis of "light cone" quantization schemes [8], [14], [15]. In terms of the momenta A 
and r associated with the null coordinates u and v the energy-momentum relation 
is At = n = ^- from (5.2). If we use the standard operator replacements for these 
momenta 



Then the the integral operators T UjV can be written symbolically in terms of 
the momentum operators as follows, 



T u = j, r„ = ^, (6.2) 



The expression (4.10)giving the solution can then be written as 



1G 



(f)(u,v) 




/(«) + e 





(6.3) 



The operators ^- and ^- 



are equal to £ and A by (5.2)and play the role 



of Hamiltonian operators relative to v and u; equation (6.3)thus represents the 
operation of a pair of Hamiltonians evolving a solution from the initial values on 
the pair of null surfaces [11], [12], [13], [14], [15]. 

Note from the above discussion that initial value data on both null surfaces are 
necessary for the evolution of a solution. From the discussion of the exponential 
solutions it is also clear that the constant term in the generating functions plays 
an important role in the evolution. 

7. Acknowledgements 

We would like to thank Prof. Stanley Brodsky for his continuing support as 
well as thank the SLAC Theory Group for its hospitality. 



17 



REFERENCES 

1. R. Courant and D. Hilbert, Methods of Mathematical Physics, Vol.2 (In- 
terscience, John Wiley and Sons, New York, 1962); see chapters 5 and 
6. 

2. R. J. Adler, M. Bazin, and M. M. Schiffer, Introduction to General Relativity, 
2nd ed. (McGraw Hill, New York, 1975); see chapter 4 for the characteristics 
of the Maxwell equations and chapter 8 for those of the vacuum Einstein 
equations; the Kerr solution and the black hole null surface are disucssed in 
chapter 7. 

3. R. Penrose, Gen. Rel. and Grav. 12, no. 3, 225 (1980). 

4. H. M. Zum Hagen and H.-J. Seifert, Gen. Rel. and Grav. 8, no. 4, 259 (1977). 

5. R. P. Kerr, Phys. Rev. Letters, 11, 237 (1963); the Kerr solution representing 
a rotating black hole is built from a metric containing a null vector; see also 
ref. 2, chapter 7. 

6. S. Hawking, Commun. Math. Phys. 43, 199 (1975). 

7. N. D. Birrell and P. C. W. Davies, Quantum Fields in Curved Space (Cam- 
bridge University Press, New York, 1982); see chapter 8. 

8. S. J. Brodsky, G. McCartor, H. C. Pauli, and S. S. Pinsky, Particle World 
3, no. 3, 109 (1993). 

9. O. C. Jacob, Phys. Lett. B347, 101, 1995. 

10. O. C. Jacob, Phys. Rev. D50, 5289 (1994). 

11. O. C. Jacob, Phys. Rev. D51, 3017 (1995). 

12. C. Itzykson and J.-B. Zuber, Quantum Field Theory (McGraw Hill, New 
York, 1980); see chap. 1. 

18 



13. C. Morosi and L. Pizzocchero, J. Math. Phys. 35, 2397 (1994). 

14. V. O. Sovoliev, J. Math. Phys. 34, 5747 (1994). 

15. J. Losa and J. Vives, J. Math. Phys. 35, 2856 (1994). 



19 



